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Abstract 

Some exactly solvable potentials in the position dependent mass background are generated whose 
bound states are given in terms of Laguerre- or Jacobi-type X\ exceptional orthogonal polynomials. 
These potentials are shown to be shape invariant and isospectral to the potentials whose bound state 
solutions involve classical Laguerre or Jacobi polynomials. 
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1 Introduction 

In recent years, quantum mechanical systems with a position dependent mass have attracted a lot of 
interest due to their relevance in describing the physics of many microstructures of current interest, 
such as compositionally graded crystals [I], quantum dots [2], 3 He clusters [3], metal clusters [I] 
etc. The concept of position dependent mass comes from the effective mass approximation [H]-[5] 
which is an useful tool for studying the motion of carrier electrons in pure crystals and also for 
the virtual-crystal approximation in the treatment of homogeneous alloys (where the actual poten- 
tial is approximated by a periodic potential) as well as in graded mixed semiconductors (where the 
potential is not periodic). The attention to the effective mass approach stems from the extraordi- 
nary development in crystallographic growth techniques which allow the production of non uniform 
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semiconductor specimen with abrupt heterojunctions. In these mesoscopic materials, the effective 
mass of the charge carriers are position dependent. Consequently the study of the effective mass 
Schrodinger equation becomes relevant for deeper understanding of the non-trivial quantum effects 
observed on these nanostructures. The position dependent (effective) mass is also used in the con- 
struction of pseudo-potentials which have a significant computational advantage in quantum Monte 
Carlo method [HI [10]. This large variety of applications explains the growing interest in constructing 
solvable cases of effective mass Schrodinger equation. Recently, several of such cases were obtained 
by means of point canonical transformation (PCT), Lie algebraic techniques and super symmetric 
quantum mechanical methods \T\ - 32 . 

In this letter, our objective is to construct new exactly solvable potentials in a position depen- 
dent effective mass background whose bound state wavefunctions can be written in terms of recently 
introduced [33J [M] Laguerre or Jacobi type X% exceptional orthogonal polynomials. (It must be 
mentioned here that similar study has been done by Quesne [35] in the constant mass scenario) 
Subsequently, we shall show that these exactly solvable potentials are shape invariant by using su- 
persymmetric quantum mechanics (SUSYQM) technique in position dependent mass case. SUSYQM 
allows us to obtain isospectral partner potentials of these exactly solvable potentials. The motivation 
for doing this, apart from enlarging the class of exactly solvable potentials for position dependent 
mass Schrodinger equation, comes from the fact that in different areas of possible applications of 
low dimensional structures e.g. quantum well, quantum dot, there is need to have energy spectrum 
which is predetermined. Specifically, in the quantum well profile optimization isospectral potentials 
are generated through supersymmetric quantum mechanics. These are necessary because a partic- 
ular effect such as intersubband optical transitions in a quantum well, may be grossly enhanced 
by achieving the resonance conditions i.e. appropriate spacings between the most relevant states 
and also by tailoring the wave functions so that the (combinations of) matrix elements relevant for 
this particular effect are maximized J3B]. This is particularly important for higher order nonlinear 
processes. Optimization of simple stepwise constant profiled quantum wells has been considered 
quite some time ago [37J, while the optimization of continuously graded structure requires more 
sophisticated techniques like supersymmetric quantum mechanics and inverse spectral theory [38J. 

The organization of this article is as follows. In section 2. we have generated two new exactly 
solvable potentials for one dimensional position dependent mass Schrodinger equation whose wave- 
functions involve Laguerre or Jacobi-type X\ extended polynomials via point canonical transforma- 
tion approach(PCT) [39, 40J. In section 3. we have shown that the new potentials are shape invariant 
in the context of supersymmetric quantum mechanics in the position dependent mass background. 
Section 4. is kept for summary and discussions. 
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2 Generation of new potentials via PCT 



The general position dependent mass (time independent) Schrodinger equation, initially proposed 
by von Roos [H] in terms of three ambiguity parameters r, s, t such that r + s + 1 = — 1 is given by 

where the effective potential 

s + 1 M" M' 2 
Veff = V(x) + -±- w - [r(r + s + l) + s + l] 1 ^ (2) 

depends on some mass term M(x). Here prime denotes differentiation with respect to x. M{x) is 
the dimensionless form of the mass function m(x) = moM(x) and we have taken h = 2rriQ = 1. 
Let us find the solution of ([1]) of the form 

^{x)=f{x)F{g(x)) (3) 

where f(x),g(x) are two function of x to be determined and F(g) satisfies second order differential 
equation of the type 

fl + Q{g) ^ + R{g)F ^ (4) 



Now in order to be physically acceptable, ip(x) has to satisfy the following two conditions : 
(i) Like quantum systems with constant mass it will be square integrable over domain of definition 
D of M(x) and ip(x) i.e., 

\il) n {x)\ 2 dx < oo 



(ii) The Hermiticity of the Hamiltonian in the Hilbert space spanned by the eigenfunctions of the 
potential V(x) is ensured by the following extra condition [25] 

hMx)| 2 ;Q 

^M{x) ^ 

at the end points of the interval where V(x) and ip n {x) are defined. This condition imposes an 
additional restriction whenever the mass function M(x) vanishes at any one or both the end points 
of V. 

Substituting ip(x) from equation ([3|) in equation (P) leads to 

d?F fg" 2f> M'\dF ( f M M'f \ 

Comparing equations (jU) and ([5]) we get 

Q(9(*)) = ^ + 7F-W 



(6) 
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The first of Eqns. ([6]) gives 



/(x)oc r~ exp i i f 9(x) Q( t ) dt ) ( 7 ) 




while the latter leads to the equation 

F _ v = jT 3 (9"Y , 9' 2 ( E ldQ Q 2 \ M" 3M* 

e// 2Mg> AM\g')^M\ 2 dg A J 2M 2 ^ 4M 3 {) 

For equation (|8|) to be satisfied, we need to find some functions M(x), g(x) ensuring the presence of 
a constant term on the right hand side of equation ([8]) to compensate E on its left hand side and 
giving rise to an effective potential V e ff(x) with well behaved wavefunctions. 

Now in PCT approach there are many options for choosing M(x), for example M(x) = Xg' 2 (x) 

[TU129], M = Xg'(x) , M = — — [42] , A being a constant. Here we choose M(x) = Xg'(x) so that 

9 \ x ) 

equation flSJ) reduces to 

*-^-$(-5f-?) 



Example 1 

st we cons 

j. For this polynomial 



First we consider F n (g) cx where n = 1, 2, 3, a > is the Laguerre type X\ polynomials 



9(9 + a) g g + a 

For these values of Q{g) and R{g) equation (9) becomes 

F _ v = j ( 2an + a 2 -a + 2 \ g' ( a 2 - 1 \ 2g' g' g' 

eff Xg\ 2a J g 2 \ 4A J X(g + a) 2 Xa(g + a) 4A 1 j 

Taking = C, where C is a constant, a constant term can be created on the right-hand side of 
the above equation which will correspond to E on the left-hand side. But C must be restricted to 
positive values in order to get increasing energy eigenvalues for successive n values. The solution 
of the above-mentioned first order differential equation for g(x) leading to a positive mass function 
reads 

g (x) = e~ bx , M(x) = e~ bx , - oo < x < oo (11) 

where CX = —b, b > 0. This exponentially behaved mass function is often used in the study of confined 
energy states for carriers in semiconductor quantum well \27\ I28j . This mass function was also 
found to be useful for understanding transport properties through semiconductor heterostructures 
(produced by the recent crystallographic growth techniques) which is indispensable for the prediction 
of the performances of these samples . In this context, this mass function has been used to compute 
transmission probabilities for scattering in abrupt heterostructures |46j which may be useful in the 



4 



design of semiconductor devices [37]. 
Now equations (f7|) . (PTOl) . and (flTl) give 



E = h — | I/; + - + 2n - 2 j -!- \\ 

a 



eff 



bx i „— fez 



(a 2 - l)e bx + e 



+ 



+ 



4 Va(l + ae fe:r ) (1 + ae 1 



+ ^0 



where Vq is a constant. Thus we obtain energy eigenvalues and eigenfunctions as 



Em = bf[m + 



a + 1 



+ - + V 

a 



ip m (x) =7V„ 



exp[-\{{a + l)bx + e bx }] f(a) , 



m = 0,1,2... 



a + e bx 

where we have reset quantum number n = m + 1 and M m is the normalization constant given by 



(12) 

(13) 
(14) 



Mr, 



b ml 



[m + a + l)r(m + a) 



1/2 



It is to be noted that the square integrability condition(i) stated earlier does not impose any additional 
restriction on a but for condition(ii) to be satisfied, a should be greater than — |. In figure 1 we 
have plotted the mass function M(x) given in (jlip . potential V e ff given in (|12p and square of first 
two bound state wavefunctions. 




^ M (x) 




\ \ 




\ 









Figure 1: Plot of the potential V e ff (solid line) given in equation ([12]) . square of its first two bound 
state wave functions ^(x)! 2 (dashed line) and | -01 C^) 1 2 (dotted line), for the mass function M(x) 
(long dashed line) given in equation (|11|) . We have considered here b = l,a = 2. 
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Example 2 : 

Next we choose F n (g) oc P^^\x), where Pn*'® , n 
type X\ polynomial. For this polynomial, we have 



1, 2, 3, a, (3 > — 1, a ^ [3 is the Jacobi 



Q(g) 
R{g) 



2W-a) 



(a+P+2)g-(P-a) 

1=P (f3-a)g-(P+a) 
(P-a)g-(n-l)(n+a+/3) _ Q3-a) 2 

T^g 2 (p-a)g-(P+a) 



(15) 



For these Q(g) and R(g) equation dH) becomes 

E ~ Veff 

where 



' v A ig + A 2 A 3 g + A A 



Ac 



i-g 2 



+ 



A fi 



A, 



_ /? 2 -a 2 
(3 2 -a 2 



(l-<? 2 ) 2 (j3-a)g-(J3 + a ) [((3 - a)g - (J3 + a)] 2 

/3 2 + a 2 
2a/3 



(16) 



1 



A 2 = n 2 + (/? + a - l)n + -{(/? + a) 2 - 2(/3 + a) - 4} + 



Ai 



/3 2 + a 2 



2a/3 



To generate a constant term on the right-hand side, we suppose — - 



constant. Consequently we obtain 
g(x) = tanh(ax) , 



A(l-9 2 ) 



Ci > 0, where C\ is a 



M(x) = sech 2 (ax) 



OO < X < oo. 



(17) 



where C±X = a, a > 0. This asymptotically vanishing mass function depicts a solitonic profile |48j . 
Recently, this mass profile has been used in position dependent mass Hamiltonians of Zhu-Kroener 
[44j and BenDaniel-Duke [43] type and interesting connection was shown [48] between the discrete 
eigenvalues of such Hamiltonians and the stationary 1-soliton and 2-soliton solutions of the Korteweg- 
de-Vries (kdV) equation that match with the mass function up to a constant of proportionality. 
Also, in dealing with position dependent mass models controlled by a sech 2 -m&ss profile, it was 
demonstrated [45] that in the framework of a first order intertwining relationship, such a mass 
environment generates an infinite sequence of bound states for the conventional free-particle problem. 
Now using ([7j). (|16p . (|17p and setting quantum number n = m + 1, we obtain the new potential , 
energy eigenvalues and corresponding bound state wavefunctions as 



v; 



eff 



((a 2 - l)e 2ax + ((3 2 - l)e~ 2ax ) 



+ 



4(a-/3)(a-3/?) _ 8/9(a - f3) 
a((3 + ae 2ax ) a{(3 + ae 



2 1 a + (3 
a m H 



a + P + 2\ 2 fj3 a 2 + (3 2 -2 , 
m+^— ]+a (- i )+V 



(19) 



0+1 



1p m (x) = M„ 



'l-tanh(ax)) * (1 + tanh(ax)) pM> {tanKax)) 



a + j3 + (a — P)tanh(ax) 



where Vq is a constant and the normalization constant M m is given by [33 



AC 



a(a - (3) 2 m\(2m + a + (3 + l)r(m + a + /? + 1) 
2 a +P- 1 (m + a + l)(m + (3 + l)r(m + a)r(m + (3) 



m = 0,1,2, 



1/2 



(20) 
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An additional restriction a,(3 > — 1/2 is to be imposed to satisfy condition(ii) stated before whereas 
the square integrability condition does not require any extra restriction on the parameters. In figure 
2 we have plotted the mass function given in (|17p and potential V e ff given in (|18j) and square of its 
first two bound state wavefunctions. 
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Figure 2: Plot of the potential V e ff (solid line) given in equation (|18p. square of its first two bound 
state wave functions ^(aOI 2 (dashed line) and | , 0i(j;)| 2 (dotted line), for the mass function M(x) 
(long dashed line) given in equation (fTT|) . We have considered here a = .2, a = 2,/3 = 2.5. 

It should be mentioned here that for the suitable choices of the constants Vo and Vo, the ex- 
pressions within the square bracket of the obtained potentials given in (|12p and (|18p . coincides 
with previously obtained potentials in [32] using classical Laguerre or Jacobi polynomial. For the 
same choices of Vo and Vo, the potentials obtained here also become isospectral with the previously 
obtained potentials in [42] . 



3 Supersymmetric quantum mechanics approach 

Defining 

Aib = _L_^ + Bib and A*ib = ~4~ I -t= I + Bib 
y/M dx dx VVM/ 

where B{x) = — jfjtr is the superpotential, the Hamiltonian of equation ([!]), can be factorized in 



/m i>o 
the following way 

H eff = AU = -— (—^—) — + V eff (21) 
" dx \M{x) J dx n K ' 

It's supersymmetric partner Hamiltonian is given by 

H^-AA* = _^(_L_)| + Fw/ (22) 
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where V e ff and Vi te ff are the supersymmetric partner potentials given by 



v, 



eff 



B 

Jm 



+ B 1 



Vl, ef f = V ef f + 



2B' 



(23) 



These two potentials will be called shape invariant if they satisfy the condition |18l [49] , 

V 1>e ff(x, ai) = V eff (x, a 2 ) + R(ai) 



(24) 



where a\ is a set of parameters , a 2 is some function of a\ and R(a\) is independent of x. In case 
of unbroken supersymmetry , the energy spectrum and wavefunctions of two such shape invariant 
effective mass potentials are related by 



£<>//) =0 , E£*ff>=E$P ,n = 0,1,2,.. 



and 



(l,e//) 



'W.'^^n±i_ ; ^ff = "/" ,n = 0,l,2... 



E 



(eff) 
n+l 



p(l,e//) 



(25) 
(26) 



Now the superpotential -B(x) for the potentials obtained in the Example 1 and Example 2, can 
be written as 



/ % ^ r / \ — — — 

B(x) = - (a + \)e 2 — e 2 



6 e2 



a(a + l)e 2bx + (2a + l)e 6x + 1 



and 



B(x) = — [(a — 0)cosh(ax) + (a + + 2)sinh(ax)] 



+ 



2a(a - 0) 



(27) 



(28) 



[(a + 0)cosh(ax) + (a — 0)sinh{ax)} [(a + + 2)cosh{ax) + (a — /3)sm/i(aa;)] 

respectively. 

Now using ([23]) and (f2T|) we obtain 



(a 2 - l)e bx + e~ bx + 



„bx 



a(l + ae bx ) (1 + ae 



bx\2 



h 2 r 



a(a + 2)e b:c + e~ bx + 



Jjx 



+ 



(a + 1)(1 + (a+ l)e fex ) (1 + (a + l)e te ) 2 



L 




- + 






a. 


Si 


^ a 







(29) 



a + 1 



Also using ([23]) and ([28]) we obtain 



(a A - l)e 2ax + (0* - l)e-' zax + 



2aa; , 4(a-/3)(a-3/3) 80(a - 0) 



2 i 



^ ^2a/3 + 2a + 2/3 + 2 + — I 



4/A 
a / 



a(0 + ae 2ax ) a(/3 + ae 2ax ) 2 



(30) 



(31) 



a(a + 2))e 2a * + (/?(/? + 2))e~^ + 



2(KC , 4(a-/?)(a-3/3-2) 



8(/3 + l)(a-/3) 2 



4 V(a + l)((a + l)e 2aa; +^ + l) 



(a + l)((a + l)e 2aa; +/?+l) 
a + 1 



(32) 
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The potentials (V e ff) obtained in ([29j) and (f3Tj) are same with the potentials (fl~2|) and (fT8|) for 
Vo = -6 2 (^ + ~) and Vo = ~t ( 2 «/ 3 + 2a + 2 ^ + 2 + ¥) respectively. 

Prom the equation (|29|) and (|30p we observe that the the potential V e ff and it's super symmetric 
partner potentials V\^ e ff satisfy the following relation 

V heff (x,a) = V eff (x,a + l)+b 2 (33) 

Also from equation (f3Tj) and ([32]) it is clear that the potentials V e // and Vi, e // are related by 

Fi je// 0, a, /?) = y e //(x, a + 1, /? + 1) + a 2 (a + /3 + 2) (34) 

From the above two relations we observe that the potentials obtained in Example 1 and Example 2 
satisfy the condition (|24l) . So we conclude that the new potentials are shape invariant. 
Now using (68), (73) of [33] and ([26]) we have derived the eigenstates of the partner potential of the 
potential (fl~2]) , as 

WW,) oc >(.->*) , m = ,l,2... (35) 

(a + l)e 2 + e 2 

and using (42), (47) of [34J and (f26j) we obtain the eigenstates of the partner potential of the potential 
CEU) as 

l+j3 1+ci 

^ff)(x) oc ( 1 + tanh ( ax )) 2 (l-tanh(ax)) 2 p(«+W+i) (tanh(ax)) m - 1 2 

^ m iXj ^ (a - flsinhiax) -(a + (3 + 2)cosh(ax) m+1 ^ an ^ ax >> > m - U, 1, 2... 

(36) 



4 Summary and Outlook 

We have obtained exactly solvable potentials for position dependent (effective) mass Schrodinger 
equation whose bound state solutions are given in terms of Laguerre or Jacobi type X\ exceptional 
orthogonal polynomials. As mentioned earlier, the obtained potentials are the generalizations (by 
some rational functions) of the previously obtained potentials |42j whose bound state solutions in- 
volve classical Laguerre or Jacobi orthogonal polynomials. The method discussed here can be used 
for other choices of the function g{x) in order to generate other type of exactly and quasi exactly 
solvable potentials for the one dimensional Schrodinger equation with position dependent mass. We 
have shown that these potentials are shape invariant and are isospectral to the previously obtained 
potentials in position dependent mass background whose solutions are given in terms of classical 
Laguerre or Jacobi type orthogional polynomials. Though the origin of such isospectrality in the 
constant mass scenario has recently been shown in |50] , it will be worthwhile to study the origin of 
such isospectrality in the position dependent mass background. As to the possible physical applica- 
tions of our obtained potentials in the position dependent mass background, let us mention that it 
will be interesting to use supersymmetric quantum mechanics to generate isospectral potentials that 
depend on a specified number of scalar parameters by multiple deletion and restoring of some levels 
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of the original potential in the position dependent mass scenario so as to make them suitable for mul- 
tiparameter optimization of optical nonlinearities in semiconductor quantum wells. This particular 
problem was dealt in ref |51j . taking the constant effective mass m*. 
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